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Abstract : The present paper provides a multi-period contagion model in the credit risk field. 
Our model is an extension of Davis and Lo's infectious default model. We consider an economy of n 
firms which may default directly or may be infected by other defaulting firms (a domino effect being 
also possible). The spontaneous default without external influence and the infections are described by 
■ not necessarily independent Bernoulli-type random variables. Moreover, several contaminations could 
i be required to infect another firm. In this paper we compute the probability distribution function of 
the total number of defaults in a dependency context. We also give a simple recursive algorithm to 
D ' compute this distribution in an exchangeability context. Numerical applications illustrate the impact 
of exchangeability among direct defaults and among contaminations, on different indicators calculated 
from the law of the total number of defaults. We then examine the calibration of the model on iTraxx 
data before and during the crisis. The dynamic feature together with the contagion effect seem to 
have a significant impact on the model performance, especially during the recent distressed period. 



Keywords : credit risk, contagion model, dependent defaults, default distribution, exchangeabil- 
ity, CDO tranches 

1 Introduction 

(N ' 

I> ' The recent financial crisis marked the need for paying more attention to the systemic risk which can 
. partially be the result of dependence on many factors to a global economic environment. A tractable 
\^ and common way of modeling dependence among default events is to rely on the conditional indepen- 
dence assumption. Conditionally on the evolution of some business cycle or macroeconomic related 
■ factors, defaults are assumed to be ind eper ident. However, as shown by s ome empirical studies such as 
O I Das. Duffie. Kapadia and Saita ( 200?! )^ or lAzizDOur and Gieseckel (ffl), the latter assumption seems 
^ I to be rejected when tested on historical default data. An additional source of dependence, namely 
the chain contagion effect, is observed and requires the construction of contagion models which would 
be able to explain the "domino effects": a defaulting firm causes the default of another firm which 
^ I infects another one etc. Additional en ipirical eyidenc e s and economic analyses of the default conta- 
gion e ffect can be found in for example iBoissav (|2006l l. |jorion and Zhan3 (|2007l l or lJorion and Zhana 



To capture these effects in a realistic way multi-period contagion models have to be considered. 



In this paper we present an extension of Davis and Lol 's model, presented in 2001 in their paper 



"Infectious Defaults" by the authors. They propose a new way of modeling the dependency between 
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defaults through "infection". The original idea applied to "CBO" {Collateralized Bond Obligation) is 
that any bond may default either directly or may be infected by any defaulting bond. Even though 
in their paper, the authors model contagion effects in a bond portfolio, their study may be easily 
applied to any elements subject to a contagion risk (for example firms, persons, buildings, etc.). In 
our paper, the application will focus on credit risk portfolio modeling. In the original model, the direct 
defaults are described by Bernoulli-independent and identically distributed random variables denoted 
and the infections are also independent and identically distributed random variables denoted Y^^ , 
all BernoulH distributed {i G N, j G N). 



The present paper provides a model more general than Davis and Lo's. We propose a multi-period 
model and we further relax the assumption of independence and identical distribution of the direct 
defaults and infections. This fact will allow the insertion of possible additional dependencies and to 
take into account the assets diversity in the portfolio. It will be also possible, for example, to consider 
different direct default probabilities or default dependencies related to the presence of a common 
economic factor. 

Moreover, compared to Davis and Lo's model in which only directly defaulting bonds can infect 
others, our model allows taking into account the domino effect which can exist between the bonds, firms 
etc. Thus in the model presented here, the firms can default because of a chain reaction, phenomena 
which is often a reason for financial crises. 

Finally, the infection way will be generalized here in order to consider the situations where several 
infections are necessary to generate a new default : it is possible for example to consider the case where 
only the default of several firms of a sector (and not only one default) can "destabilize" other firms. 
These situations can thus limit the contagion phenomena, but also can (together with the contagion 
probabilities) allow the introduction of possible critical mass of defaults predisposed to make the crisis 
more sudden (for example if a large number of infections is necessary to generate other defaults and 
if the infection probabilities are large). 

The model can be used for example for epidemic studies, accidents in chain, or to introduce in an 
original manner a dependence structure which could translate the impact of these contaminations. 



Co ntagion models were introduced to the credit risk field bv lDavis and Ld (j200lh and lJarrow and Yu 
(|200lh . Since 2001, the financial contagion has become a problem which drew the attention of 
many authors. JioO^) proposed an extension of Jarrow and Yu's model. Davis and Lo's in- 
fectious model was also studied by many au thors who proposed different mode ls having as starting 
point Davis and Lo's model (see for exarn ple Egloff. Leippold and Vanini ( 200?! ) . Rosch. Winterfeldt 



( 20081 ). Sakata. Hisakado and Morj ( 2007)). All along years, the contagio n phenomena was modeled 



using Bernoulli random variables (IDavis and Ld (l200lh). copula fun ctions (jSchonbucher and Schubert 
2001)), interacting particle systems ( Giesecke and Weberl ( 2004)), incomplete iii f ormation models 



( Frev and Runggaldiei ( 2008 ) ) or Mar kov chains (see for example Schonbucher ( 20061 ). Graziano and Rogers 
tmd ) or iKraft and SteffensenI tmii )). As far as the risk management of synthetic CDO tranches 
is concerned, Ma r kov c h ain contagion models h ave also b een i nvestigated by several papers such 




goncerneu, Ma r Kov c n am contagion models n ave also b een i nvestigated b y several papers suci 
Van der Voort ( 2006 ). Herbertsson and Rootze n (2 0o"^'). Herb ertsson ( 20071'). Frev and Bac khaul 



) , iLopatin and MisirpashaevI (j2007l ) , lArnsdorf and HalperinI (j2007l ) , ICont and Mincal (120081'). ICont. Deeuest a 



Frev and BackhausI (l2008l\.lDe Koch. Kraft and SteffensenI (l2007l 'l iEpple. Morgan and Schloeg. 



ihaud 
:hloegll 



2009) amo ng others. The hedging issue for CDQ t ranches is also addressed bv lLaurent. Cousin and Fermanian 
2003) and Cousin. Jeanblanc and Laurent ( 20091 ) in the class of Markovian contagion models. As a 
first step among possible applications of the Davis and Lo's multi-period model, we illustrate the 
model tractability in terms of the pricing of synthetic CDO tranches. Since the computation of CDO 
tranche spreads involves the computation of expected tranche losses at several time horizons, such a 
study could not have been done using the original static version of the model. We then discuss the 
model performance in terms of fitting liquid CDO tranche spreads. 
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In the paper Infectious Defaults, Mark Davis and Violet Lo present a model where the bonds 
may default either directly or may be infected by any defaulting bond (see Davis and Lo . 200 ll ). The 



authors consider n bonds, n G N*. Let be the set of possible indices of these bonds, = {1, ..,n}. 
For a bond i, i G 0, the random variable X* is equal to 1 if the bond defaults directly, otherwise. 
Even if = 0, bond i can be defaulted by infection (contagion). A default infection takes place if a 
Bernoulli random variable = 1. Thus the bond i is defaulted (directly or by infection) if = 1, 
with 

= X' + {1 - X')'ff\ 



In lDavis and Lol 's paper, the infection is modeled in the following way : 



= 1 - Y[ [1- X^Y^' 

where Y^^, i,j £ Q are Bernoulli random variables. Consequently, at the end of the period, a bond i is 
infected if at least one X^Y^^ = 1, j E 0, j ^ i, i.e. at least another bond j defaults directly {X^ = 1) 
and infects the first bond {Y^'^ = 1). 

In the following, will denote the binomial coefficient(0 ^ /c ^ n) : = "l^); ■ 

Davis and Lo stated and proved the following theorem : 



Theorem 1.1 ( Davis and Lol 's result) 

ien 

where = /(I -p)"^'=(l - g)'=("-'=) + ^ Clp^(l - p)"^^(l - (1 - g)*)*^-^(l - ?)*("-'=). 



fc-i 



i=l 

In this paper, we propose to relax the following assumptions of the original model : 

• Al. The random variables {X*, i G il} are independent and identically distributed and V i G ^2, 
F[X' = 1] =p,pG]0,l[, 

• A2. The random variables {Y^^ , G 0,^} are independent and identically distributed and 
yi,jen,j^i,F[Y'^ = l] =q,qe]0,l[, 

• A3. The infection occurs if at least one X^Y^^ = 1, j G ^2, j ^ i. 

• A4. An infected bond can not infect another bond (one period model). 



The outline of the present paper is as follows : we present the model and we fix the assumptions 
under which we will state our theorems (Section [2]). We present some life insurance tools used for the 
characterization of the defaults' number law (Section 13. ip . In the following Section [3.21 we give the 
defaults' number law in a one period model. Our result allows finding, under the random variables 
independence assumption, Davis and Lo's result. In Section 13. 3[ we generalize the model of the 
previous section and we present a multi-period model. Finally we give some numerical applications 
(see Section [3D including an illustration of parameters effects on the loss distribution and a discussion 
on the ability of the model to fit iTraxx CDO tranche quotes. 

2 The Model 

We consider an economy of n firms (n G N*). Let be the set of possible indices of these firms : 
= {1, ...,n}. We observe the firms along a given time interval divided into T periods. When the 
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firm i (i E ri) defaults directly during the period t (i € {1, ...,T}), then the random variable XI is 
equal to 1, otherwise X'l = 0. Even if = 0, the firm i can be defaulted during the period t : it may 
be infected by another defaulting firm. We use a binary random variable to denote this infection. 
This variable is equal to 1 when the firm i does not default directly but is infected by other firms 
during the period t, = otherwise. Let {Z^, i € 0} be a sequence of discrete random variables 
with possible values {0, 1}. For any i Q, Zl = 1 if the firm i is defaulted (directly or by infection) 
at the end of period t, Zl = otherwise. In our model, the random variables {Zl, t = 1, ...,r} are 
obtained by the following recursive relation : 

Zl = zi_, + (1 - zi_,)[xi + (1 - xim, 2<t<T,ien, 

Z{ = X{^{\-X\)'€\. 

Hence Z\ = 1 if the firm i has been declared in default at the end of period t — 1 {Z\_y = 1) or if 
it defaults directly {XI = 1) or by contagion {'^l = 1) during the period t. 
We introduce the following notations : 

Notation 2.1 For every t G {1, . . . T}, we denote by : 

• Qt the set of the firms declared in default up to period t : Qt = {i ^ — ^} > 

• the set of the firms which did not default in the previous t periods : Ft = {i G Vt, Z\ = 0} = 

• the number of spontaneous defaults without external influence occurred at time t : 

• Nt the cumulated number of defaults occurred up to time t : Nf = J2ieO. ■ 

For every t G {1, . . . T}, the random variable Y^^ ^ ((j, i) G VL^) with possible values {0, 1} represents 
the infection during the period t of the firm i by a defaulting firm j. Consequently the variables ^\ 
(i G 0) have the form 

where / : {0, ..,n} — )• {0,1} and Ft represents the set of the defaulting firms susceptible to infect 
the firm i. Here we suppose that card{Ft) = g{Nt-i, N^), where 5 : — >■ {0, ..,n}. For example, 
in the original model Fi = (thus only the firms which default directly may infect others) and 
f{i) = Ij^i i.e. at least a contamination causes the default. In our model, we can imagine the case 
where the firm i defaults only if it is infected by a given number of firms, for example f{i) = li^2 (i-e. 
at least two infections are necessary to generate a new default). Moreover, the firms which can infect 
the firm i during a given period t, are those which were directly defaulted during the same period t, 
or those declared in default (directly or by infection) before this period. 

The aim of this paper is to study the law of Nt under the following assumptions : 

Assumption 1 (Temporal independence of direct defaults) 

The random vectors Xt = {Xt^---iXt)j t ^ {!)••• j^^}; o,f^ mutually independent, but a dependency 
exists between the components of each vector. 

Assumption 2 (Temporal independence of exchangeable infections) 

The random vectors Yt = (Y^^^, ^/^) •••) ^"")i t G {1, . . . , T}, are mutually independent and for any t, 
the random variables (j, i) G fi^j are exchangeable, independent of {XI, t = 1, . . . ,T, i G $1}. 

Under these assumptions, we suppose that all joint distributions of (X^ , . . . , X^"), t € {1, . . . ,T} 

and of {Y^^^\ . . . , Y^^"^ ^) are known, where Y^^^ is a permutation of the set Y-l'' , {i,j) G O^. 
We will sometimes use the following assumption which is a particular case of Assumption [T] : 
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Assumption 3 (Temporal independence of exchangeable direct defaults) 

The random vectors Xt = i^ti ■■■j^t')! i G {1, ■ ■ ■ ,T}, are mutually independent, but the components 
of each vector are exchangeable random variables. 

In the particular case of infinite exchangeability assumption, thanks to De Finetti's Theorem, this 
assumption is perfectly adapted to the situation where the direct defaults are dependent because of a 
"hidden" common factor (for example the state of the economy). 

At last we will sometimes consider the particular case of independent and identically distributed 
(i.i.d.) direct defaults and i.i.d. contaminations : 

Assumption 4 (i.i.d. direct defaults and i.i.d. contaminations) 

The random variables {X'j:, i G fi, t £ {1, . . . ,T}} are independent, all Bernoulli distributed with pa- 
rameter p and {i,j) G f^^, t G {1, . . . ,T}| are also independent random variables, all Bernoulli 
distributed with parameter q. 

Compared to Davis and Lo's infectious model, we propose a multi-period model where {XI, i G 0} 
and {i,j) € are no more independent and identically distributed random variables. More- 

over, we do not consider just the particular case where f{i) = Ij^i and the firms which can infect are 
not necessary those which default directly during the same period. Thus, we take into account the 
"domino effect", i.e. the firm i defaults, it infects the firm j which infects the firm k etc. 

To generalize the mode of infection, one can for example consider, f{i) = lj^2 instead of f{i) = 
Ij^i. In other words, two contagions are required to cause an indirect default. This feature of our 
model allows to weaken the effect of contagion. Such a possibility can be useful for large portfolios. 
Indeed, it is important to generalize the mode of infection. Let us imagine that the number of 
firms increases. Then, if all the parameters of the model remain identical, too many infections could 
occur. Other recent generalizations were proposed and they introduce different solutions to generalize 
the mode of infection in order to better represent the credit risk markets reality (see for example 
Sakata. Hisakado and Mori . 200?! ). Furthermore, a generalization may allow to better model the 



critical thresholds up to which the infections start to have a real effect. Different modes of contagion 
will be compared in numerical applications (see Subsection 14. 2|) . 

3 Results 

The aim of this section is to present the main results of the paper together with some tools used to 
prove them. 

3.1 Life Insurance Tools 

In this section we present some useful tools to characterize the defaults number distribution. 

We introduce the coefficients of order k for the set {XI, i G F} where F C il, card{r) > k, k N 

Definition 3.1 Let F C il. For all t G {1, . . . , T}, the coefficient of order k (k < card{T) ) for the set 
{XI, I G F}, denoted iik,ti^), is defined as 

Mfe,i(r) = Yl P [^1" = 1 n ... n = l] , l<k< card{T), 

^card{r) ji<j2<-<jk 

/^o, i(r) = 1 (including j/F = 0). 
The ^ symbol means to sum all C^^rdir) possible choices of k different elements taken among 

jl<j2<--<jk 

the elements of the set F. 
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Remark 3.1 // the random variables {X^, i £ are exchangeable, then for allT C ^ we have 



xl = lf^...r\X^ 



l<k< card{T). 



In the particular case where {XI, i S 0} are independent random variables, all Bernoulli distributed 
with parameter p, then fif^^t = p^. 

A classical problem in the actuarial field is to compute the number of survivors after a given period, 
among a group of different persons with independent and identically distributed survival rates. We 
will see that these tools can considerably simplify the studied credit risk models. 

Lemma 3.1 (Waring formula) Let Xl,...,XJ^ be n dependent Bernoulli random variables and let 
T be a subset of {1, ... ,n} such that card{T) = m < n. Then 



m—k 



ik<mCl ci_,(-i)V,+fe,t(r). 

3=0 



This result was proved bv iFellerl (|l968l ) (chapter IV. 3, pa ge 106). An ex tension of this formula is 
Schuette-Nesbitt Formula often used in actuarial science (see iGerberj (jl995l ). chapter 8.6, page 89). 



Remark 3.2 Let X} , X^ be n exchangeable Bernoulli random variables and let Mx = {i : XI = 1] . 
Then for all C {1, . . . ,n} such that card{E^^) = k we have the following equalities : 



.1=1 



1, x^^^ 



XI 



C^P [Mx = Ek]. 



Since the variables jl^/"', [i,]) S f^^j are exchangeable, then conditionally to [Nt-i^N^Y the 

variables i G fij are also exchangeable. We introduce the coefficient of order k for these two sets 

of random variables. 

Definition 3.2 For all T dVi^ and for all t G {!,... T}, the coefficient of order k (k < card{T) ) for 
the set \y^\ (i,j) € pj, denoted Xk,t, is defined as 



Xk,t =P 
Ao,t=l. 



y: 



a(l) 



in...ny; 



a(k) 



1 



k>l, 



In the same way we define 

a,t(9(n, /)) =P [^j = 1 n ... n = i\Nt-i = u, = I 

CoAdiuJ)) =1 (including if g{u,l) = 0). 



A; > 1, 



Every variable is a function of j € -Ftj which are exchangeable random variables. The 



<7(1) 



law of {'to'I, . . . , "^J*) depends only on card{Ft) and on the law of {Y^ 
establish a relationship between the coefficients Xk^t and ^,k,t{z), z <n. 

Proposition 3.2 (Joint lav^^ of ("TJ . . . '^^)) 
The joint law of {^\ . . . is given by 



Y' 



a-{card{Ft)) 



). We can thus 



^k,tiz) 



/(5:y/^) = in...n/(^yf) = i 

j<^Ft j£Ft 



card[Ft 



, 0</c<n, {)<z<n-k. 
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1. In the particular case where f{i) = li^i, we have 



k zi 

a,t(^)=5:ci5:c^,(-ir+"A„,,. 

2. For any function / : {0, . . . , n} ^ {0, 1}, we have 

zk — 'y 

7=0 i=0 

where r]k,z{l) = E /(7i)C2'%_i_^(7 - 71), r/i,^(x) = lx<zf{x)C% and rio^x) = 
7ie{o,...,2} 

71 <7 

This proposition is proved in Appendix. 

Remark 3.3 // {y/^ (i, j) £0?, t£ {1, . . . ,T}] are independent random variables, all Bernoulli 
distributed with parameter q, then X^^t = CLnd 

• in the particular case where f{i) = (,k,t{z) = (1 — (1 — qY)^ , 

zk 

• for any function / : {0, . . . , n} ^ {0, 1}, ik,t{z) = ^ Vk,z{'y)l"' " qf'''"'. 

7=0 

3.2 One Period Model 

The aim of this section is to characterize the defaults' number distribution in the case where T = 1. 
For the sake of clarity, all subscripts 1 are omitte d, so that = ^"1, X* = X|, = 

It is a one period model which generalizes the Davis and Lol 's results when the direct defaults and 
the infections are not necessary ind ependent. In t he particular case where T = 1, the random variables 
{Z*, i ^ are obtained as in the Davis and Lol 's model : 

= x' + {i- x')'€\ i e 

where = /(^ card{Ft) = g{N^) (with g{x) = g{0,x)). Let us recall that in the original 

model, the authors consider only the particular case where f{i) = li>i. 

Under Assumption [U the random variables {X^, i G 17} are no more independent. We suppose 
that all the joint laws of (X^, . . . , X"") and ^ , . . ^ y^C"")) are known. 

One of the main results of this paper is the following theorem which gives the defaults' number 
law. 

Theorem 3.3 (Defaults' number law when the defaults are not i.i.d. - One period model) 

Under Assumptions [I\ and\^ ifT = l, then the defaults' number law is given by : 

r n—r n—k 

p [AT = r] = c; ^ cZ-A-^Tic+r-Kiim) E Ctfc(-l)V.+fc(^^)- 

fc=0 a=0 i=0 

Proof: Using the law of alternatives, we get 



V[N = r] = Yv\N = r\N'^ = k 



k=0 



N^ = k 
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On the one hand, by LemmaEH P = k = ^ C^_;.(-l)Vj+fc(^^)- 

j=0 

On the other hand 



N = r\N 



D 



k 



+ {1- X')'^' 



.1=1 



N'' = k 



= r-k = k 



N = r\N 



D 



C^ZlEaZoC^l-ri-ina+r-kAgik))- 



□ 



where A = {i e n : X' = 0} . 

Since card{A) = n — k, Lemma l3 . 1 1 gives P 

The equahty C^C^l^ = C^C^ concludes the proof. 

We easily deduce the law of N when {X*, i € f^} is a set of exchangeable variables. 

Corollary 3.4 (Defaults' number law when the defaults are exchangeable - One period model) 

Under Assumptions \^ and O assume that the variables {X^, i £ 0,] are exchangeable. Then the 
defaults' number law is given by : 

r n~r n—k 

F[N = r] = ciYl E c^^-ri-mc+r-kAgm E cL.(-i)V,+fc, 

k=0 a=0 j=0 



where the coefficients fij are introduced in Remark \3. 

As a consequence of Theorem 13.31 Remarks 13.11 and 13.31 is the following result which gives the 
defaults' number law in the particular case where {X^, i € ^2} are independent random variables, all 
Bernoulli distributed with parameter p and {Y^^ , € 0,^} are also independent random variables, 

all Bernoulli distributed with parameter q. 

Corollary 3.5 (Defaults' number law when the defaults are independent - One period model) 

Under Assumption\^ assume that {X^, i G fi} are independent random variables, all Bernoulli dis- 
tributed with parameter p and {Y^^ , {i,j) G il^} are also independent random variables, all Bernoulli 
distributed with parameter q. Then the defaults' number law is given by 



P[iV 



+r—k 



k=0 



a=0 



where ^u{z) = E Vu,zi'y)q'^ {I 

7=0 



\2'U — 7 



Under the assumptions of Corollary 13.51 s-lso assume that at least one infection causes a new 
default (i.e. f{i) = Then, using Remark 13.31 and taking g the identity function we find the 

Davis and Lo's result. 

3.3 Multi-period Model 

Generalizing the model presented in the previous section, we consider now that the time interval is 
divided into several periods. In this model the default indicator is described by the relation : 



zi = zu + (1 - + (1 - 2<t<T,ien, 

Z\ = X{ + (\-X{)'€\. 

We are interested in computing the law of the defaults' number at the end of every period t, 
t £ {!,..., T}, T > 1. Another main result of this paper is the following theorem : 



(1) 
(2) 
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Theorem 3.6 (Defaults' number law when the defaults are not i.i.d. - Multi-period model) 



Under Assumptions {1\ andl^ the defaults' number law is given by : 

F[Nt = r]= P[®t = 9t], 

etcn 

card{9t)=r 

where Qt is the set of firms declared in default up to t (see Notations \2.1\) and where 

r 

M=o 0t_icet 

card{9t~i)=u 



J2 E p(^t^ ^ - - E ci_,{-iy^j+r-u-mAu, m), 

m=o Mtcet-et-1 j=o 

card{Mt)=m 



p{A, = P 



V i G A XI = 1 ety ieB X? = 



with card{9t) = r and card{9t-i) = u, u < r. 

The first term P [Qi = 6i] may be easily computed using the recursive formula for P [Qi = 9i\Qq = 
Proof: We split 

F[Nt = r]= E ^[®t = et]. 

etcn 

card{9t)=r 

By the law of alternatives 

r 

«=o et_ic6»t 

caTd{dt--i)=u 

It thus remains to compute the term P [0^ = 9t\Qt~i = ^t-i] when 6t~i C 6t C il., card{9t) = r and 
card{9t-i) = u, u < r. 



p[et = 9t\e 



t-i — pt-i 



yie9t- 9t-i, xf + (1 - xi)'^i = 1 ety i en-9t, xi = '^i = o\et-i 



Using Remark 13.21 and then Lemma EHJ we obtain that P [Qt = ^t|0i-i = ^t-i] is equal to 



en-et-i-Mt 



r — u — m\Nt-i = u, NP = m 



C 



r—u—m 
n—u—m 



E E p{Mt,n-9t-i-Mt) 
m=o Mtcet-et-1 

card{Mt)=m 
r—u n—r 

: E E pi^t, n - 9t^i - Mt) E ci_,{-iyij+r-u-mM^^^))- 



□ 

If the direct defaults are not i.i.d., the law of Nt is given by a recursive relation (see Theorem 
3.6p . However, this relation can induce an important load of calculation in numerical applications. 
That is why in the numerical applications we will consider the case where the random variables 
{XI, i G ri, t G {1, . . . ,T}} satisfy the Assumption [3] and are exchangeable. 

The following result is a particular case of Theorem 13.61 
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Theorem 3.7 (Defaults' number law when the defaults are exchangeable - Multi-period model) 

Under Assumptionsl^ and El the defaults' number law is given by 

k=0 7=0 a=0 j=0 

The proof of this theorem rests on the following lemma shown in Appendix : 

Lemma 3.8 Under Assumptions \^ and\^ for any 6t and 9t~i such that 6t^i C 9t C il., card{9t) = r, 
card[Ot~i) = k, k < r, the following relation is true : 

P [Nt = r\Nt-i = k] = C;itP [0* = St\Qt^i = <^t-i\. 
Proof of Theorem 13.71 Using the law of alternatives and then Lemma \'i.'6\ we obtain : 

r 

P [AT, = r] = ^ P [Nt = r\Nt-i = k]P [Nt^i = k] 

k=0 

= ^ p [iVi_i = k]c:;-_\F [Gi = dt\et-i = 0t-i], 

fc=0 

for any 9t and 9t-i such that 9t-i <Z 9t dVt, card{9t) = r, card{6t-i) = k, k < r. 
By Theorem I3.6[ 



r—k n—r 

P [e* = 9t\Qt-^ = 9t^^] = Y, J2 PiMtM- 0t-i - Mt) Ci_,{-iy(,+r-k-^Mk,l))- 

7=0MtC6»t-6»i_i j=0 
card{AIt)='y 

However, under Assumptions [5] and El p{Mt, Q — 9t-i — Mt) = Y^a=o~'^ Cn-fc-7(~^)"/^7+a. *■ remark 
that ^ 1 = C^_^ which concludes the proof □ 

Mtcet-dt-i 

card{Mt)='y 

We easily deduce the law of Nt when all the variables are mutually independent. 
Corollary 3.9 (Defaults' number law when the defaults are independent - Multi-period model) 

Under Assumption assume that both {Xi} and {y/^j are i.i.d. Bernoulli random variables of 
respective parameters p and q. The defaults' number law is given by 

r r—k n—r 

P [AT, = .] = ^ P [Nt., = k]c:;-_i y: cj_,p'ii-pr-''-^ Y ci_A-iy+%+r-k-,A9ik,i)), 

k=0 7=0 j=0 

zu 

where ^^A^) = J] 7?,,,(j>^(l - g)^""^'. 

j=0 

Moreover, if f{i) = li^i, then 
P [iV, = r] = ^ P [Nt^i = fc]C;it E - p)"-'-^(l - (1 - g)9(*^'^))^-'=-T(l - q)9ik,^)(--r) _ 

k=0 7=0 

Using the last equation of Corollarv 13.91 for t = 1, g{k,^) = 7, and Nq = almost surely, we find 
the Davis and Lo's result. 
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4 Numerical applications 



In this section, we provide some numerical applications of the previous extension of Davis and Lo's 
model. 



We first investigate the effect of exchangeability, among direct defaults a nd among infections, on the 
total d efau lts' number evolution. Thi s kind of impact is also dealt with by iDenuit. Dhaene and Ribas 
(l200lh and lOenuit. Lefevre and Ute"^ tooi ) in a life insurance framework, and it appears that some 
quantitative measures could be notably affected by dependencies among the considered random vari- 
ables. Here, we will be concerned about the impact of such dependencies on the first moments of total 
defaults' number Nt and on the survival function of Nt at some points. 



We then study the ability of the model to price synthetic CDO tranches. Interestingly, the first 
application proposed in the original paper by Davis and Lo focuses on the rating of Collateralized Bond 
Obligations. However, to the best of our knowledge, no calibration procedure of the Davis and Lo's 
contagion model has been tested so far on liquid tranche quotes. In this framework, direct defaults are 
considered to be Bernoulli mixtures and contagion events are assumed to be independent. Using this 
specification, we compare the calibration performance of the model on several years of iTraxx data 
before and during the crisis. 



4.1 General settings 

4.1.1 Algorithm in the exchangeable case 

Let us consider a framework where both direct defaults and infections are described by the following 
two sets of exchangeable random variables {X^, i £ Q] and jl^*"'' (^li) ^ f^^j- So that Assumptions [2] 
and [3] hold. Furthermore, suppose that the random vector laws do not change over time. As a 
consequence, in order to lighten further notations, we omit the t subscript in coefficients S,k,t{z), fik,t 
and Xk,t, which become (,k{z), fik and Afc. The presented algorithms could be easily adapted in the 
case where those quantities vary over time. 

From Theorem 13.71 we have seen that the law characterization of total defaults' number was 
requiring the computation of the set of coefficients {Ckiz), k G {0, . . . , n} , z £ {0, . . . , n — k}}. In the 
general case where / is not necessarily the indicator function f{x) = l^^i, given initialization values 
that are provided in section [XT} and given the Proposition [321 this set of coefficients can be computed 
by the following algorithm: 

Proposition 4.1 (Algorithm for {^^(-2), k,z £ N}) The coefficients {CfcC-^), A;, z G N, + z ^ n} can 

he computed by: 

• Initialization: for z varying from to n, r/o,z(0) = 1, Coiz) = 1. 

• For k varying from 1 to n, 

— For z varying from to n — k, 

* Uz) ^ 0, 

* For 7 varying from to kz, imin < — max(0, 7 — z), imax < — min(7, {k — l)z), and 

*max 

zk—'y 

ik{z) ^ a(^)+r/M(7) E Cifc_^(-l)'V7- 

j=0 
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The law of the total defaults' number, with exchangeable direct defaults and exchangeable infec- 
tions, in the multi-period model, results directly from Theorem 13.71 and is given by the following 
algorithm: 

Proposition 4.2 (Algorithm providing the law of A^^ given the coefficients Ck{z)) Once the co 
efficients {(,k{z), k = 0, . . . ,n z = 0, . . . ,n — k} calculated, the law of Nt is given by: 

• Initialization: for r varying from to n, P [Nq = r] = lr=o- 

• For t varying from 1 to T, 

— For r varying from to n, 

r r~k n—k—'y n—r 

P[Nt = r] = Y^P [iV,„i = kjCl-l E Cj-k E C^fc_,A^,+„ E Ci_A-iy^%+r-k-Mk 

k=0 7=0 a=0 j=0 

4.1.2 Beta-mixture model 

We have chosen a framework enclosing the one by Davis and Lol . where contagions in a given pe- 



riod are issued from the direct defaults during this period, so that the function g : $7 x — )• N 
is given by g{k, 7) = 7. Both direct defaults and infections are represented by sets of exchange- 
able random variables. Therefore it becomes necessary to specify the sets {Xj, j G {O, ...,n^}} 
and {nj, j € {0, . . . ,n}}, which characterize these direct defaults and infections. We assume that 
Xi, . . . , Xn are part of an infinite sequen ce of exchangeable Bernoulli random variables. Then, by De 



Finetti's Theorem (see iDe Finettil ()193ll )). there exists a random variable Gx, with values in [0, 1] 
and a cumulative distribution function Fq^, such that P X^ = 1, . . . ,X^ = 1 = f^ 9^dFQ^{6). The 
direct defaults X^ , . . . behaves as if the Bernoulli parameter p was becoming a hidden common 
random variable Qx- As a consequence, we can write = E \ qQ . We wih take E [Qx] = p m. order 



to use the same notations as the ones used by iDavis and Ld We will suppose that the parameter 



a\ = V [Qx] is given. In the same way, we will suppose that variables (i, j) € are Bernoulli 

random variables, with common hidden parameter Gy. We will write g = E [Gy] and cjy = V [Gy] 
and we suppose that these parameters are known. We will use for both Gx and Gy some laws of a 
same family, distributed on [0,1], where moments are analytically given. For j € N, coefficients 
and \j correspond therefore to: 

/ij = m{j,p,ax), (3) 
Aj = m{j,q,aY), (4) 

where the function m{j, r, a) gives the moment of order j of a random variable G from the considered 
family of laws, with mean r and standard deviation a. 

We have chosen to use Beta distributions for Qx and Gy, which are parameterized by their mean 
and standard deviation (it could have been possible to choose other laws, like for example the Ku- 
maraswamy's law). 

Moments m(j, r, cr) are easy to compute, but we need to consider two separate cases: 



When cr^ = 0, the hidden random variable G (which may represent Qx or Gy) is a Dirac mass, 
and we find the i.i.d. case for which 

m(j, r, 0) = r-* . (5) 

When fj^ > 0, for a Beta distribution with mean r = E [G] and variance o"^ = V[G], the 
parameters Or^a and Pr,cT of the distribution can be easily deduced from r and cr^, and we then 
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get all moments of order j, j E N (F is here the Euler Gamma function): 



ar,a = r ( _ 1 ) , (6) 



= (1 - r) (^^^^ - l) , (7) 

4.2 Effect of model parameters on the dynamics of loss distributions 

All numerical values that we used for the parameters in this section remain to be specified. We grouped 
them into the following Definition 14.11 

Definition 4.1 We define 4 reference models. These models have some shared characteristics: they 
all consider 10 firms (n = IQ), on a time interval divided into 10 periods (T = lOj, with a direct 
default probability p = 0.1 and an infection probability q = 0.2. The four models are distinguished on 
the nature of direct defaults or infections, i.i.d. or not (depending on the two parameters ax and uy ) 
and on the infections number that is required to cause a default, depending on the function f . The 
specificities of the four models are presented hereafter: 

• model 1: ax = 0, cy = 0, f{x) = l^^i (i.i.d. case, one required contamination). 

• model 2: ax = 0, ay = 0, f{x) = 1x^2 (i.i.d. case, two required contaminations) . 

• model 3: ax = 0.2, cry = 0.2, f{x) = l^^i (exchangeable case, one required contamination). 

• model 4-' = 0.2, cry = 0.2, f{x) = 1x^2 (exchangeable case, two required contaminations). 

Starting from the same model assumptions, we have checked that, the coefficients ^kiz), k G 
{0, . . . , n}, z € {0, . . . ,n — k} computed either from formula 1 or formula 2 of Proposition 13.21 give 
the same loss distribution P [Nt = k],k G fl when f{x) = Ix^i- In this last case, if furthermore 
ax = (Ty = and t = 1, we have checked that we get by these two formulas the same values for 
P [Nt = k],k & Q as the one computed with Davis and Lo's result. 



Figure 1: Evolution of E [A^t] (left) and V [Nt] (right) as a function of t, for all the models described in 
Definition 14. II The curves plotted with dotted lines correspond to the case where direct defaults are i.i.d. 

The evolutions of expectation and variance of Nt as a function of t, for all of the four models that 
are described in Definition [21 are illustrated in Figures [TJ We observe that the mean of the number of 
defaults distribution is an increasing function of time. This is obviously consistent with the situation 
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where more defaults are expected in a larger period of time. In the exchangeable case (models 3 and 
4), direct defaults as well as infections arise from probabilities that depend on hidden factors. If the 
distribution of Nt is intuitively more dispersed, the impact of such a dispersion on the mathematical 
expectation is not trivial. For the mean of Nt, we can observe on Figure [T] (left side) that the expec- 
tation of Nt is only slightly modified when choosing exchangeable random variables, but this change 
is very small in the case where two contagions are required to cause a default, i.e., when f{x) = 1x^2 
(models 2 and 4). This suggests that when the contagion effect is weakened as this is the case here, 
the mean of the loss distribution is mainly explained by the mean of the direct defaults, which are 
the same in models 2 and 4. When considering the variance of Nt plotted in Figure [T] (right side), we 
observe more evident results. The variance of the loss distribution is a hump-shaped function of time. 
It is quite intuitive that the dispersion level increases as time goes by since the loss process can reach 
a larger part of its support with higher probabilities. However, when the expected number of defaults 
attains a critical threshold, the number of expected surviving names then decreases, leading at some 
point in time to a decrease in the dispersion level of the loss distribution. Let us remark that this 
occurs earlier for models 1 and 3 for which the expected number of defaults is greater. Moreover, the 
presence of a variation in the hidden random variables Qx and 0y naturally implies a larger increase 
of the variance of Nt. 



Figure 2: Evolution of P [Nt ^ 6] and P [N = 10] as a function of t, for all the models described in Defini- 
tion [ITT] The curves plotted with dotted lines correspond to the case where direct defaults are i.i.d. 

Figure [2 shows the evolution of P [Nt ^ k] , for respective default number k = 6 (left side) and 
A; = 10 (right side), among the ten considered firms {n = 10). With the chosen parameters, we observe 
that the impact of exchangeability is not systematically the same. The growth of the volatility of Nt 
does not always lead to the increase of P [Nt ^ k] , especially when this last probability is larger in the 
independence case. When we consider less frequent events, as illustrated on the right side of Figure [2] 
where the probability that all firms default is plotted, we find, on the contrary, that the impact of 
exchangeability is far more explicit. This seems logical, since the increase of the dispersion level of 
Nt, due to hidden parameters @x and By, tends to emphasize the tail of the distribution. This is 
clearly a required model behavior especially when one has in mind some applications to the pricing 
of synthetic CDO tranches. Indeed, one can think of the base correlation breakdown in March 2008, 
when no implied correlation can be found for CDX super-senior tranches. 

4.3 Calibration on liquid CDO tranche quotes 

This section examines the fit of the model to tranche spreads of the 5-years iTraxx Europe main index 
at two fixed points in time, namely August 31, 2005 and March 31, 2008. These two dates have been 
chosen to facilitate comparison of model calibration before and during the credit crisis. Let us recall 
that synthetic CDOs are structured products based on an underlying portfolio of reference entities 
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subject to credit risk. It allows investors to sell protection on specific risky portions or tranches of 
the underlying credit portfolio depending on their desired risk-profile. We concentrate our numerical 
applications on the most liquid segment of the market, namely CDO tranches written on standard 
CDS indexes such as the iTraxx Europe main index. As illustrated in Figure El the iTraxx Europe 
main index contains 125 investment grade CDS, written on large European corporations. Market- 
makers of this index have also agreed to quote standard tranches on these portfolios from equity or 
first loss tranches to the most senior tranches. Each tranche is defined by its attachment point which 
is the level of subordination and its detachment point which is the maximum loss of the underlying 
portfolio that would result in a full loss of tranche notional. The first-loss 0-3% equity tranche is 
exposed to the first several defaults in the underlying portfolio. This tranche is the riskiest as there is 
no benefit of subordination but it also offers high returns if no default occurs. The junior mezzanine 
3-6% and the senior mezzanine 6-9% tranches are less immediately exposed to the portfolio defaults 
but the premium received by the protection seller is smaller than for the equity tranche. The 9-12% 
tranche is the senior tranche, while the 12-20% tranche is the low-risk super senior piece. 



iTraxx Europe 
Main 



125 CDS, 
Investment Grade 



Second Super 
Senior 



First Super 
Senior 



Senior 
Senior Mezzanine 
Junior Mezzanine 

Equity 



100% 



20% 

12% 

9% 

6% 

3% 

0% 



Figure 3: Standardized CDO tranches on iTraxx Europe Main. 



Since individual credit exposures are all the same among names in a standard CDS index, the 
key quantity involved in the pricing of standard tranches is the cumulative loss per unit of nominal 
exposure defined by 



Lt 



1=1 



1 - Ri)Zl, 



where Ri denotes the recovery rate associated with name i. We recall that referred to as the 
time-t default indicator of name i, the latter quantity being defined by Equations ([1]) and ([2]) as in the 
previously presented multi-period extension of the Davis and Lo's model. The cash-fiows associated 
with a synthetic CDO tranche with attachment point a and detachment point b (a and b in percentage) 
are driven by losses that affect the tranche, i.e.. 



L 



{a,b) 



(Lt - a)+ - {Lt - by 



We refer the reader to Cousin and LaurentI ( 20081 ) for a detailed description of cash- flows associated 
with synthetic CDO tranches. The model that we consider for the calibration of CDO tranche quotes 
belongs to the class of models presented in section 14.1.21 More specifically, we assume that, for any 
t > 0, direct defaults XI, i = 1, . . . , n are Bernoulli mixtures where the common hidden parameter is 
a Beta-distributed random variable with mean p and variance a^- For the sake of parsimony and in 
order to avoid numerical difficulties when n is large, we consider a version of the model where only 
three parameters are needed. We consider that, for any t >0, contagion events Y^'-' , 1 < i, j < n are 
independent Bernoulli random variables with the same mean q (uy = 0). We consider here the case 
where only one contagion is required to cause a default. Using the latter restrictions, the discrete- 
time contagion model is stationary and entirely described by the vector of annual scaled parameters 
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a. = {p,ax,Q)- Let us recall that the computation of CDO tranche spreads only involves the expecta- 
tion of tranche losses at several time horizons. In the case where recovery rates are the same across 
names and equal to a constant R, it is straightforward to remark that the current cumulative loss is 
merely proportional to the current number of defaults. The algorithms described in Subsection 14.1.11 
can then be used properly to compute CDO tranche spreads. 

Let us denote by sq, si, §2, S4, S5, sq the market spreads associated with (respectively) the CDS 
index, [0%-3%], [3%-6%], [6%-9%], [9%-12%] and [12%-20%] standard iTraxx tranches and by so{a), 
si{a), S2{a), 54(0), 55(0), SQ{a), the corresponding spreads generated by the contagion model using 
the vector of parameters a. The calibration process aims at finding out the optimal parameter set 
a* = (p*,(T^,g*) which minimizes the following least-square objective function 



For both data sets, in order to analyze the calibration efficiency in a deeper way, we have compared 
the global calibration with three alternative ones, where some of the available market spreads were 
excluded from the fitting. Here are the calibration procedures we have considered: 

• Calibration 1: All available market spreads are included in the fitting, 

• Calibration 2: The equity [0%-3%] tranche spread is excluded, 

• Calibration 3: Both equity [0%-3%] tranche and CDS index spreads are excluded, 

• Calibration 4: All tranche spreads are excluded except equity tranche and CDS index spreads. 

In all calibrations the interest rate is set to 3%, the payment frequency is quarterly and the recovery 
rate is i? = 40%. We provide in Table [T] and Table [2] model spreads and optimal parameters resulting 
from the four benchmark calibration processes performed on August 31^*, 2005 iTraxx quotes. Table E] 
and Table m present the same types of outcomes but for calibrations on March 31*^*, 2008 iTraxx quotes. 

As can be seen from Tables [T] and El the calibration of the three parameters on all market spreads 
is rather disappointing. This is not surprising, especially given the poor calibration performance of 
standard factor models during the crisis, when the fit is achieved on all tranches and index quotes. 
Moreover, we observed that calibration results are fairly dependent on the initial starting point spec- 
ified in the minimization routine. We do not guarantee that the presented outcomes correspond to a 
global minimum of the objective function. However, one can note that the calibration error decreases 
when we subsequently exclude the equity tranche quote and the index quote in Calibrations 2 and 3. 
Interestingly enough, the fit on 2008 data is nearly perfect in Calibration 3 (see Table ED compared 
with Calibration 3 on 2005 data (see Table [T]) . Unsurprisingly, as illustrated by results from Calibra- 
tion 4, the fit on equity and index spreads only is perfect. We have checked that this is actually the 
case for all tranches when they are jointly fitted with the index. This can be seen as a fundamental 
required behavior of the model since we try to fit three parameters on two market quotes. Let us recall 
however that the base correlation framework had some difficulties to fit the super-senior tranches in 
the same period (March 2008). As for the optimal parameters, we observe that the one- year mean 
and dispersion level of direct default events (parameters p* and cr^) are higher in 2008 calibrations 
compared with 2005 calibrations. This can be explained by the drastic shift in credit spreads between 
these two periods. However, we can note that the optimal contagion parameter q* has the same or- 
der of magnitude in both periods except for Calibration 1 in 2008 where it appeared to be insignificant. 

To summarize, in both data sets, the best calibrations were obtained when both the CDS index 
and equity tranche spreads were excluded. Moreover, all calibrations tested on the 2008 data seem to 
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0%-3% 3%-6% 6%-9% 9%-12% 12%-20% index RMSE 



Market quotes 


24 


81 


27 


15 


9 


36 




Calibration 1 


20 


114 


7 


1 


1 


29 


0.64 


Calibration 2 




62 


32 


18 


6 


8 


0.41 


Calibration 3 




55 


29 


18 


7 




0.22 


Calibration 4 


24 










36 






Table 1: iTraxx Europe main, August 31*^*, 2005. The market and model spreads (in bp) in the four 
calibrations and the corresponding root mean square errors. The [0%-3%] spread is quoted in %. All 
maturities are for five years. 





* 

P 




* 

q 


Calibration 1 


0.0016 


0.0015 


0.0626 


Calibration 2 


0.0007 


0.0133 


0.0400 


Calibration 3 


0.0001 


0.0025 


0.3044 


Calibration 4 


0.0014 


0.002 


0.1090 



Table 2: iTraxx Europe main, August 31*^*, 2005. Optimal parameters a* = {p* , cj^, g*) resulting from 
the four calibration procedures (cry = 0). 





0%-3% 


3%-6% 


6%-9% 


9%-12% 


12%-20% 


index 


RMSE 


Market quotes 


40 


480 


309 


215 


109 


123 




Calibration 1 


28 


607 


361 


228 


95 


75 


0.25 


Calibration 2 




505 


330 


228 


112 


68 


0.20 


Calibration 3 




478 


309 


215 


109 




0.002 


Calibration 4 


40 










123 






Table 3: iTraxx Europe main, March 31*^*, 2008. The market and model spreads (in bp) in the four 
calibrations and the corresponding root mean square errors. The [0%-3%] spread is quoted in %. All 
maturities are for five years. 







P* 




q* 


Calibration 


1 


0.0124 


0.0886 





Calibration 


2 


0.0056 


0.0518 


0.0400 


Calibration 


3 


0.0012 


0.012 


0.2688 


Calibration 


4 


0.0081 


0.0516 


0.0589 



Table 4: iTraxx Europe main, March 31*^*, 2008. Optimal parameters a* = {p*,ax,q*) resulting from 
the four calibration procedures (ay = 0). 
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outperform calibrations achieved on the 2005 data. Even though additional investigations have to be 
achieved, the latter result suggests that the dynamic feature together with the contagion effect seem 
to have a significant impact on pricing performance, especially during the recent distressed period. 

5 Conclusion 

In this paper, we examine a multi-period extension of the original Davis and Lo's contagion model. 
The model allows for both spontaneous or direct defaults and defaults resulting from a chain reaction 
where defaulted names at the beginning of a given period can infect other firms by contagion at the 
end of this period. Moreover, the number of defaulted names required to trigger a contagion effect 
can be specified as a model input. At last, random variables governing direct defaults and contam- 
inations may depend on common macroeconomic factors. In a rather general setting, we provide a 
recursive algorithm for the computation of the number of defaults distribution at several time horizons. 

In numerical applications, we restrain ourselves to the case where, between each time interval, both 
direct defaults and contamination events are described by mixtures of Bernoulli variables mixed with 
a Beta-distributed random factor. We first illustrate the effect of model parameters on the dynamics 
of loss distributions. In particular, we investigate the factors dispersion impact and the contagion 
effect on the variance and the quantiles of the loss distribution. We then consider the fit of the model 
parameters on iTraxx data before and during the credit crisis. This allows to exploit the dynamic 
feature of the model and illustrate its tractability when the number of reference entities is large (and 
equal to 125). In both data sets, the best calibrations were obtained when both the CDS index and 
equity tranche spreads were excluded with a perfect fit in March 2008. In that case, we can remark 
that the contagion probability is meaningful. Surprisingly, all calibrations tested on the 2008 data 
seem to outperform calibrations achieved on the 2005 data. Even though additional investigations 
have to be achieved, the latter results suggest that the dynamic feature together with the contagion 
effect seem to have a significant impact on pricing performance, especially during the recent distressed 
period. 

Eventually, this study paves the way to additional applications such as the management of coun- 
terparty risk arising in CDS or CDO tranches transactions. 



6 Appendix 



In this part, we give the proofs of Proposition 
Proof of Proposition 13.21 



and Lemma 



1. In the particular case of infinite exchangeabili ty assumpt i on, th e first part of Proposition 13.21 
is a consequence of De Finetti's Theorem (see De Finetti ( 193ll )). Indeed, since the variables 
,j) G fi^l are exchangeables, there exists a random variable Gy e [0, 1], such that, 

conditionally to Gy, random variables ^Y^-',[i,j) G are i.i.d., Bernoulli distributed with 
parameter Gy. As a consequence 



6,t(^) =E 



j€Ft 
card{Ft)=z 



E 



(1-(1-Gi 



\z\k 



'1-G> 
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Noting that E [(1 - Gy)'^] = E 



j&A 
card(A)=zi 



we conclude using Lemma [XT! 



In the case of finite exchangeabihty assumption, De Finetti's Theorem may not be apphed. 
That's why it is necessary to use another reasoning : 

Let m > 1. Let us define for k € {1, . . . , m — 1}: 



d1(^) = p 



^ y/' > 1 n ... ^ > in ^ yi.— fc+i = q n • • • ^ y/'"^ = o 

j&Ft j£Ft j&Ft jGFt 



y: = ... = e = 



card{Ft) = z 



It follows that for all k {1 < k < j) : 



card{Ft) = z 



Since the variables |y/"', (i, j) € il^j are exchangeables, we see that for /c G {1, . . . ,m}, 



m=A; 



By induction on /c = 1, . . . , m, we get: 

=1 - ckg(^) + cieg(^) + ... + {-\fG\itl^z). 

Since the variables {yI^ ■, {hj) G are exchangeables, then for all r (1 < r < /c) , we conclude 
using Lemma 13.11 : 



E^. 



E 

Q = 



E C2j(-l)"Aa,t . 



2. Let f-^ = {ien:i< z, f{i) = 1}. Then 



iuM = E P 



^y/^ = 7in...n^y/'= = 7. 



card{Ft) = z 



Furthermore random variables Iv^"**, (i, ^) G are exchangeable. So, using Remark 



■ ■ - ^z p 

/-i7lH h7fc 

-1 ^fcz 



Ei;'^^^'^ = 7i + --- + 7. 



7iv>7fcG/; 

fcz (-(71 /~j7fe 

E E /(7i).../(7.) ^ - p 

7=0 71 H l-7fe=7 

7i>---,7feG{0,.--,z} 



p7 



From Lemma l3.ll 



kz 



7 



fcz— 7 

cL E cL_,(-i)^A,+.,. 
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As a consequence 

7=0 j=0 

where i^uAl) = E /(Ti) • • • /(7/c)Cr ...CI". 

71 H l-7fe=7 

7iv,7fce{0,...,2} 

A simple calculation leads us to the following reccurence relationship : 



%,^(7) = E /(7i)CIi%-i,^(7-7i; 



7iG{0,...,2} 

71 <7 



Proof of Lemma 13.81 From the definition of conditional probability 

Y>[Nt = r,Nt-i = k] 



□ 



V[Nt = r\Nt^i = k] 



P [Nt-i = k] 



Under Assumptions [2] and [3l the random variables {Zl, iEil,tG{l,...,T}} are exchangeable. 
On the one hand, using two times Remark 13.2^ we establish that P [Nt = r,Nt-i = k] is equal to 



_(~\k (-\r—k 
-^n^n-k 

■^k r~\r—k 



Nt = r,ZU = ■■■ = ZU = = ••• = ZU = 

Nt = r,Zl = -.. = Z'^ = ZU = ■■■ = ZU = 1,^'-+/ = ••• = ZU = 

/7I ryr ryl ryk -1 lyT + l tyn ryk + 1 ryn r\ 

— ■ • • — — — • • • — ^t-l — J-, — • • • — — _ . . . _ Zj^_^ — u 



On the other hand, using the same remark, P [Nt-i = k] = C^P [Qt~i = which concludes the 

proof. □ 
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